A Green-function approach to transport phenomena in quantum pumps. 
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I. INTRODUCTION 

The impresive development in the technology of fabri- 
cation of small circuits, enabled the investigation of single 
electron transport induced by time-periodic fields [lj . In 
a recent experiment by Switkes et al Q charge transport 
through a quantum dot is induced by means of two time- 
periodic potentials with a phase lag. This has renewed 
the interest in the study of charge pumping in open quan- 
tum systems, motivating theoretical 0, El IE El II BIH 
and experimental activity 0]. Closel y re lated phenom- 
ena are the photovoltaic effect [H Hi 

m mcso- 

scopic j unctions and charge driving in molecular ratchets 
[lj, . The main feature characterizing these effects is 
the generation of a net current as a response to a time 
dependent external field without a net static bias. 

Several theoretical treatments on quantum pumps rely 
on adiabatic approximations 0- LU- 0- LJ| ■ relevant for very 
slow potential modulations. A good amount of work 
on time-de pen dent transport is based on Floquet the- 
ory [SSHHIHEIIZI, the matrixapproach I1I1II3 
or on the transfer matrix technique [ljj. An alternative 
framework to investigate transport phenomena in meso- 
scopic devices and nanostructures is the description of 
the device in terms of tight-binding Hamiltonians and 
the solution of the problem with non-equilibrium Green 
functions. Since the proposals of Refs. [2I this 
kind of approach became widely extended in the study 
of electronic transport through a nanometric or meso- 
scopic sample, as a response to a static bias. One of the 
reasons for the success of this scheme is the fact that it 
is suited to deal with arbitrary high bias, finite tempera- 
ture, arbitrary strength of dissipation and that it can be 
extended to include many-body interactions at least per- 
tubatively. Another appealing feature of this approach 
is the possibility of combining it with the so-called "ab- 
initio" methods to describe details of the contacts and 
molecular bridges of the devices j2^| . In the context of 
transport problems with time-dependent fields, there are 
also basic proposals of this type of strategies 
but, in comparison, not so many recent developments. 
Some examples are studies on ac-driven quantum dots 
and superlattices [H El El |28| , the latter restricted to 
weakly coupled quantum wells, studies on the dynami- 



cal Franz-Keldysh effect [29| , superconducting point con- 
tacts with a time-dependent voltage 30], and conduct- 
ing rin gs t hreaded by a time-dependent magnetic flux 
EH EH- In these problems, the time-dependent 

par t of the Hamiltonian is restricted to a single point 
[Hill, bond [UHl Him or contacts [HUll, while 
approximations are introduced to deal with more gereral 
situations [H El El El 

In this work, we present a general treatment based 
on non-equilibrium Green functions to study transport 
phenomena in systems described by tight-binding mod- 
els in contact with particle reservoirs with several time- 
periodic local potentials. We derive exact equations of 
motions and present results on the transport properties 
in the special simple cases of one-dimensional systems 
with one and two time-dependent potentials. In the first 
case, some analytical expressions are available. In the 
second one, which is relevant for the experimental con- 
figuration of Ref. 0, we compute the Green functions 
numerically. 

The traditional and intuitive way to think about sta- 
tionary transport through a mesoscopic device placed be- 
tween two electrodes at different chemical potentials is in 
terms of the behavior of the density of states of the cen- 
tral system and of its environment. In a time-dependent 
problem the density of states depends on time and its 
convolution with a Fermi function does not directly cor- 
respond to the notion of occupied energy states. In our 
study we analyze the connection between the transport 
behavior of the pumps, the density of states of the en- 
vironment and the non-equilibrium spectral densities at 
the positions where the time-dependent potentials are 
applied. 

In stationary transport like that resulting as the re- 
sponse to a static bias, the carriers responsible for the 
transport process are those injected from the electrodes 
with energies between the two different chemical poten- 
tials. A remarkable property of the pumping mechanism 
is that not only electrons with energies close to the Fermi 
energy of the reservoirs contribute to the net electronic 
current. Instead, all the electrons contribute to the net 
flow. This point has been previously addressed in Ref. 
|l9j | for the problem of a time-dependent harmonic po- 
tential in an asymmetric structure. We present here fur- 
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FIG. 1: (Color online) Scheme of the system described by 
Hamiltonian Q with L = 4 reservoirs and M = 3 pumping 
potentials. The grid represents the sites of the central sys- 
tem and the boxes with stripes represent the reservoirs with 
chemical potentials fi a , with a = 1, ... ,4. 



ther details on this behavior which combines effects like 
photon-assisted tunneling, quantum interference and dis- 
sipation, sometimes giving rise to patterns that resemble 
a turbulent motion of electrons through the device. 

The control of the direction of the net current is cen- 
tral for eventual technological applications of the pump- 
ing effect. However, the complex nature of the electronic 
motion generated in a quantum pump makes the predic- 
tion of this property from a priori considerations a very 
difficult task. We discuss some operational conditions 
where the sign of the induced current can be understood 
and we also identify some mechanisms causing sign rever- 
sals. The paper is organized as follows. The theoretical 
treatment is presented in section II. Examples and results 
are presented in section III and IV. Finally, section V is 
devoted to summary and conclusions. 



where the fermionic operators a,k a denote degrees of free- 
dom corresponding to the reservoir a, which is described 
by the Hamiltonian H a . The contact of the reservoir a 
with the central region is represented by the hopping ele- 
ment w a between the reservoir and the sites l a placed at 
the boundaries of the lattice of the central system. The 
model for the latter piece is a tight-binding Hamiltonian 
with TV lattice positions, hopping elements Wij between 
pairs of sites and local energies £/. We also consider 
time-dependent potentials acting locally at M sites of the 
central system denoted by ii , oscillating with amplitudes 
Vi, frequencies Q; = nflo (being n an integer number) 
and phases Si. A scheme of the setup is shown in Fig. ^ 

In what follows, we present a closed set of equations 
to calculate the Green functions corresponding to spacial 
coordinates of the central system. As usual |35j , we define 
retarded and lesser Green functions 



G£ m (M0 = -;e(i-t0<{c*(*),4(*')}>, 

Gfc.Ct.O = i(cl(t') Cl (t)). (2) 



In the non-equilibrium formalism, Dyson equation for 
the Green function has a matricial structure which re- 
sults in a coupled set of integro-differcntial equations for 
the lesser and retarded components. Following the stan- 
dard procedure polEllEl iBl l2l| we consider ideal non- 
interacting leads for the reservoirs. The corresponding 
degrees of freedom are integrated out, defining the self- 
energies 
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II. THEORETICAL TREATMENT 

A. General model and Green functions 

We assume a general device, which consists in L reser- 
voirs and a central system driven by several time-periodic 
local potentials. We do not include the effect of many- 
body interactions in our treatment. The full system is 
described by the Hamiltonian 



L L 

H = Y Ha - Y Wa Y °L c ^ - Y w u4 c j + 

a=l a=l k a ,l a (lj) 

N M 

Y s i4 c i +Y Vicos ( nit +s ^ c l c ii +h.c, (i) 
i=i 1=1 



which depend on the equilibrium Green functions 
9k 'k — °f the free reservoirs (i.e. isolated from the 
central system). The latter can be expressed in terms of 
densities of states p a (u)) through 

k J -00 Z7r 

/CO 7 
— e- <W ^f a (cj) Pa (u), (4) 
-co 27r 

being / Q (w) = l/(e^( w_AlQ ' + 1) the Fermi function. 
We work in units where h = 1. The equations for the 
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retarded and lesser components read 

8 N 

1 = 1 

-Y, / dt 1 G^ >l (t,t 1 )E R n (t 1 -t') = 6(t-t')S m , n , 
i=i Jt ' 

8 N 
{ ~ l 8V ~ £n ~ ^')}G<,„(M') -£<?< ,(M'K„ 

i=i 

-£[/ dt 1 G^(t,t 1 )E< n (t 1 -t') 
1=1 J -°° 

+ [ dt 1 G<,(t,t 1 )tf n (t 1 -t')] = 0, (5) 

J — OO 

being V n (t) = S n>il Vj cos(Qi* + Z = 1, ...,M. 

We now present a very convenient way to calculate the 
Green functions. We start from the following integrated 
form of eqs. 

N ft ft' 

G< >n (t,t') = E / dti / dt 2 GtSM) 

iJ=l J ~°° J -c° 

x£<(i 1 -t 2 )G£ n (t a ,t'), 

G* n (t,i') = G^, n (t -0+5] / dhG^h) 

i=i Jt ' 

xV^Gl^h-t'). (6) 

The advanced Green function is related to the retarded 
one through Gf i7l (t,t') = [G^j(t' ,t)}* . The equation for 
the lesser Green function is valid for long enough t,t', 
such that no memory on the initial condition is preserved. 
The "unperturbed" retarded Green function G° mn (t-t') 
corresponds to the solution of the problem of the central 
system without time-dependent potentials coupled to the 
reservoirs. 

We now define the Fourier transform 

G% tn (t,u>)= [* dt'GKjtXy^^'K (7) 

J — oo 

with r\ = + , which leads to a linear set for the retarded 
Green function. 

G^(t, W ) = G^ n H+53|e^+^) 

i 

xg« , (t, w - n,)G? » + £ ^ e -^+ n ") 
xG^.ct.w + no^nM. (8) 

In most of the cases, this set must be solved numer- 
ically, by discretizing u> introducing high and low fre- 
quency cut-offs ±K£l , respectively. The parameter K 
depends on Qq and must satisfy that |iff2o| i s much larger 



than the absolute value of the highest frequency for which 
G^ n (u>) has a finite spectral weight. In the numeri- 
cal procedure, it must be checked that the solution of 
G^ n (t,uj) does not depend on K. 

Since G^ ln (t,u>) is a periodic function of t with pe- 
riod To = 27r/57o, it is sometimes useful to work with the 
expansion: 

oo 

„(*,")= E <W^K fcnot , (9) 

k= — oo 

being 

m , n (M) = - r dte-* kn ° f GK n (t,u;). (10) 
T o Jo 

We stress that the advantage of writing the Dyson 
equation in the form © and working with the Fourier 
transform J7J is that the exact retarded Green function 
can be evaluated from a set of linear equations irrespec- 
tively the amount of time-periodic potentials. In addi- 
tion, the present formulation in terms of eq. JSJ is very 
convenient to perform systematic expansions in powers 
of VI. 

B. Pumped current in a one-dimensional system 
coupled to two reservoirs 

We now consider the case where the pumped system 
is a tight-binding chain with hopping elements between 
nearest neighbors, which is placed between left and right 
reservoirs. The Hamiltonian reduces to 

H = H L +H R + H c (t)-w L (a\ j c l +H.c.) 

-w R (a^ R c N + H.c), (11) 

with Hc{t) denoting the Hamiltonian for the central 
piece. 

The current from the reservoirs to the central region 
can be written as 

J a (t) = 2eu> Q Re[G£,(i,t)], (12) 

with a — L and i = 1 (a = R and i — N) for the 
current flowing from the left (right) reservoir. The dc 
components of the above currents are 

1 f T ° 

J d a c = - dtJ a (t), (13) 
T o Jo 

and due to the continuity condition, they satisfy J| c = 
— J R C which allows us the write the dc current J flowing 
through the central device as J = (Jf c — J R )/2. 
Following j23| we write 

1 f T ° f* 

J=- dt dt 1 Re{\w L \ 2 [G? 1 (t,t 1 )g<(t 1 -t) + 

T JO J-oo 

Gl^t^gtih-t)] - \w R \ 2 [G*, N {tM)9i{h-t) 

+G<^(t,t 1 )g£(t 1 -t)},}, (14) 
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Alternatively, it is also possible to calculate J from the 
dc component of the current flowing trough an arbitrary 
bond (I, I + 1} of the central tight-binding chain; 

J u+1 (t) = 2ew M+1 Re[G< +1 (M)]- (15) 

Due to the continuity property, the dc component of l|15|) 
is independent of I and coincides with the result obtained 
from {nj. 

We consider zero temperature and the same chemical 
potential \i for the two reservoirs, hence / Q (w) = f{u) — 
0(/i — lo). An interesting representation of J is found 
expressing it in terms of a transmission function T{u>): 

/oo 
duof{uj)T{uj). (16) 

Evaluating Gf l+1 (t,t) from © the explicit equation for 
the transmission function reads 

T ( W ) = [ T \t{\w L \ 2 p L (cj)lm[Gf :i (t,u;)Gt +1 ^,t)} 
+ \w R \ 2 p R (u)lm[G* N (t, u)G^ tl+1 ( U , *)]}. (17) 

The above representation of J exhibits a very important 
difference between the current generated in a quantum 
pump and the stationary transport caused by reservoirs 
at different potentials. In the latter situation, the current 
is expressed as J = e J dtuT '(u>)[f — fn(u)}. Namely, 
the spectral contribution of T(lu) to the net current corre- 
sponds to states with energies between the two different 
chemical potentials. Instead, in the pump all the states 
bellow the Fermi energy of the reservoirs contribute to 
the net current. Another important difference between 
these two kinds of transport mechanisms is the origin 
of the direction of the current. In stationary transport, 
T(w) is a possitive-defined function, which is interpreted 
as the probability of tunneling, while the sign of J is de- 
termined by the bias through /l(^) — fai^)- Instead, in 
time- dependent transport without static bias, the trans- 
mission function (|17|) can be either possitive or negative 
and can also change sign as a function of uj. This function 
can be interpreted as the difference between the proba- 
bility of tunneling from left to right and the probability 
of tunneling from right to left. 

The next sections are devoted to evaluate explicitly 
T{ui) and J for the particular cases of one and two time- 
dependent potentials and to analyze in detail their be- 
havior. 



III. ONE HARMONICALLY 
TIME-DEPENDENT POTENTIAL. 

A. General considerations. 

The treatment exposed in the previous section simpli- 
fies considerably for the case of only one harmonically 



time-dependent potential (M = 1). The Hamiltonian for 
the central system reads: 

H c (t) = [Vcos(n t) - ei]c\ci. (18) 

The reservoirs L and R are placed at the left and the 
right, respectively, of the pumping center, as explained in 
the previous subsection. Our aim is the calculation of the 
net current. Using the expansion J!|J), the dc components 
of the currents flowing from the reservoirs towards the 
central site read 

/OO 7 
— /(w){-2Im[gi,i(0, w)]p a (w) 

- £)/>a(w - kn )T{Lo)\g 1A (k, lu)\ 2 }, (19) 

k 

being T(u) = \wl\ 2 Pl(u) + |wr| 2 /9r(w). It is clear that 
a necessary condition for a non- vanishing J^ c is that the 
two terms of (|19|) do not cancel one another. 

Further insight is gained by using the representation 
of the current in terms of the transmission function. In 
the present case, it is possible to solve recursively the 
set (jHJ) to calculate the retarded Green function. This 
procedure is summarized in Appendix A. After replacing 
the expressions (|A5|) in l|14|) . with N = 1, it is found 

T(u,) = -\w L \ 2 \w R \ 2 V \Gi,x{k,u>)\ 2 

TT — ' 

k— — oo 

x{p l (lo)pb{uj - kilo) - pr(w)pl{u - kVL )}, (20) 

where it becomes clear that J = when (i) the sys- 
tem is symmetric under spacial inversion centered at the 
point where the pumping potential is applied, such that 
Pl(uS) = pr(lv), and (ii) the environment of the pumping 
point can be described by flat and approximately con- 
stant densities of states p a (uj). This expression for T(lu) 
resembles the mechanism of photon-assisted tunneling. 
The resulting equation for the current has a similar form 
as that obtained within the framework of the scattering 
matrix formalism, identifying the square of the scatter- 
ing matrix in the Floquet formalism with the function 
\Gi,i(k, uj)\ 2 pl(li->)pr(uj — kflo) 0. The current can be 
viewed as the result of processes where electrons leave 
the reservoir a with probability p a {u>), interact with the 
pumping center loosing or gaining k energy quanta Qq 
with probability ex \Gi,i(k, w)| 2 , and exit to the opposite 
reservoir (3 with probability pp(u) — fcfio). The sign of the 
net current is completely determined by the structure of 
the functions p a {ui). Note that T(u>) may change sign as 
a function of w. This means that electrons with different 
energies can flow in different directions and it is the sum 
of all these contributions what determines the sign of the 
net current. 

For small Oo and V, only a few modes k contribute. It 
is natural to associate such a situation with the idea of 
adiabatic pumping. More precisely, the concept of adia- 
batic pumping applies to the regime where the character- 
istic time scale for an electron to travel across the pump 
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(proportional to the inverse of the width of the spectral 
peaks of G(k, u>) as a function of u>) is much smaller than 
fl . Adopting that definition, we see that it is possible in 
this case to have a finite J even in the adiabatic regime. 
The key is a high hybridization wr,wl of the pumping 
center with the environment, in order to allow for wide 
peaks in Q{k,Lo), meaning a short life of the electrons 
at the pumping center. For low pumping frequencies, 
(|20|l can be expanded in powers of fig and it is found 
T(lo) oc fig. 

A final interesting remark is that J oc |u>l| 2 |u^| 2 . This 
kind of behavior has already been found in molecular 
ratchets pumped by a laser field [l^l , which are modeled 
on the basis of tight-binding Hamiltonians with asym- 
metric energy profiles and sincronic pumping centers. 



B. Example 

To illustrate the discussion of the previous subsection 
we show some results of local pumping in a symmetric 
double barrier structure with the pumping potential act- 
ing at one of the barriers. The "unperturbed" structure 
has spacial inversion symmetry with respect to the cen- 
ter, which is broken by the effect of the time-dependent 
voltage. 

A scheme of the device is shown in Fig. [21 For sake of 
clarity we write down the model Hamiltonian: 



H = Hi eads + H cont -w ^ cjci+i + Hc + 

l=-N L 

1 

^2 £ i4°i + Vcos(n Q t)c\c 1 , (21) 

l=-N L 

where ei defines the profile corresponding to the barri- 
ers. We consider a two-barrier structure of height Eh'. 
E-n l = El = Eb and Ei = 0, 1 ^ —Nl, 1. We denote with 
Hieads the Hamiltonians of two semi-infinite chains which 
behave as macroscopic reservoirs and represent two ex- 
ternal leads connected to the central device. These parts 
are pictorially represented by boxes with stripes in the 
scheme of Fig. [21 and we describe them by semicircular 
densities of states with bandwidth W, 



P°{oj) = 4^1-uj 2 /W 2 Q(W - oj). 



(22) 



The term H cont describes the hopping between the semi- 
infinite leads and the central structure and has the form 
of the last two terms of the Hamiltonian , with hoping 
parameter wq. 

In the notation of Eqs. (|18|l and the right reser- 
voir corresponds to the right lead, while the left reser- 
voir corresponds to the left lead plus the double barrier 
structure with the exception of the point 1 = 1, where the 
pumping voltage is applied. The ensuing degrees of free- 
dom can be easily integrated out defining the following 




FIG. 2: (Color online) Solid red line: the function t 1 (ui) = 
\G(l, cu)\ 2 for V = 0.1, n = 0.1. Dashed black lines: The 
densities of states Pl{uj) and Pr(ui) of the left and right reser- 
voirs, respectively. A denotes the energy difference between 
two consecutive energy levels in the double barrier structure. 
A scheme of the device is indicated in the upper left corner. 



Other parameters are Eb = 1, |wo| 



0.1, W = 4, N L 
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retarded Green functions for the left and right reservoirs: 



g*(u J )=g (cj) 



duj' p°(uj') 



2tt u> — us' 



(23) 



with rj — + , and 
g&u) = 



UJ - E2 — 



CO ■ 



w - e Nr - w 2 3°(w) 

being the densities of states of the reservoirs p a (u)) = 
— 2Im[<7„ (w)]. The hoppings between the pumping center 
and the left and right reservoirs are wl = w and wr = 
wo, respectively. 

Fig. [21 illustrates the behavior of the function tk{uj) = 
|^i i i(A;,iu;)| 2 , related to the probability for an incoming 
electron with energy to to loose k energy quanta Qq for 
selected parameters, along with the densities of states of 
the left and right reservoirs. All the energies are written 
in units of the hopping parameter w, which is set to w — 
1. All currents are expressed in units where e = 1. The 
results shown in the figure correspond to a weak pumping 
amplitude (V = 0.1), where only the contribution k = 1 
is sizable. For larger V, higher modes come also into 
play. 

The resulting structure of the transmission function 
T(u>) is shown in Fig. [31for a low frequency Qq and differ- 
ent amplitudes V. The important feature to note is that 
for high u> and large pumping amplitude, this function 
can experiment several changes of sign. This indicates 
that electrons with different energies may flow along dif- 
ferent directions. The direction of the total direct current 
being the cumulative sum of all these contributions. 
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FIG. 3: (Color online) The transmission function T(ui) for 
flo = 0.1 and V = 0.1, 0.5, 1 (top to bottom). Other parame- 
ters are as in Fig. [5] 
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peaks separated in ~ fio (see bottom panel of Fig. [3Jl. 
This is translated in peaks at the edges of the plateaus 
shown in J as a function of fi as shown in the bottom 
panel of Fig. Details of the structure related to f2o is 
shown in the inset. 

For yU < 2, the chemical potential lies below the highest 
resonance of /?l(w) and the behavior of T(lo) is consis- 
tent with a preference in the flow right to left, while for 
higher [i, the flow can take place in the opposite direc- 
tion. Details depend on the frequency and amplitude of 
the pump, as well as on the degree of coupling between 
the double-barrier structure and the macroscopic reser- 
voirs. 

The behavior of the sign of J for /i < 2 is roughly 
the one expected from an intuitive adiabatic description. 
The picture that emerges is as follows: during the part 
of the pumping cycle, where the potential decreases, the 
tunneling from the right macroscopic reservoir into the 
quantum well is favored, while during the remaining part 
of the cycle, the total potential at the right barrier in- 
creases and the electrons accumulated in the well are 
pushed from the right to the left. Remarkably, within 
this range of /x, the sign of J inferred from such a simple 
description remains unchanged for high and strong V, 
where the adiabatic picture is not expected to be valid. 
The situation for higher fi is much less clear and strongly 
depends on the particular values of the pumping param- 
eters. 

To finalize, we want to remark that, due to the sym- 
metry of the problem, the situation where the pumping 
potential is applied at the left barrier, would result in 
exactly the same behavior of J, but with the opposite 
sign. As discussed in the previous subsection, a vanish- 
ing J would be obtained if the pumping potential were 
applied exactly at the center of the well between the two 
barriers. 



FIG. 4: (Color online) The dc current J as a function of 
the chemical potential of the electrodes fi. Top and bot- 
tom panels corresponds to pumping voltages V = 0.2 and 
V = 1, respectively. Plots in thin solid black, dashed red, 
dot-dashed blue and thick solid magenta lines correspond to 
fio = 0.01,0.05,0.1,0.2, respectively. The spacing between 
two consecutive levels is indicated with A. The inset shows 
details for V = 2 and Qo = 0.1. Other parameters are as in 

Fig. m 



Some results on the behavior of the net current J as a 
function of the chemical potential fi of the reservoirs are 
shown in Fig. 01 The structure of jumps and plateaus 
observed in the figure follows the patterns of resonances 
related to the peaks of t^{oS) and pl{u). The energy 
interval between two consecutive jumps in J is roughly 
the difference of energy between two energy levels A. For 
strong pumping amplitude, many k modes come into play 
in the transport process. The transmission function de- 
velops a richer structure where the spectral weight asso- 
ciated to the free electronic levels split into several side 



IV. TWO HARMONICALLY 
TIME-DEPENDENT POTENTIALS IN A 
DOUBLE BARRIER STRUCTURE. 

A. General considerations 

We now go back to Hamiltonian (|llfl and consider the 
following Hamiltonian for the central region: 

N 

H c (t) = V cos(n t + S)c\ Cl +J2eicjci 

i=i 

N 

—w y^(cjcj + i + H.c) + V cos(r2oi)cjvCAr, (25) 
1=1 

with the profile E\ — En = Eb, £; = 0, 1 = 2, . . . , N — 1, 
defining a double barrier structure. This arrangement is 
similar to the one of the experimental setup of Ref. 0] , 
where two ac potentials with a phase-lag are applied at 
the walls confining a quantum dot. 
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In the interesting case of reservoirs with wide bands, it 
is possible to find an explicit relationship between the 
current and the local spectral functions at the points 
where the pumping potentials are applied. The wide 
band limit corresponds to approximately constant densi- 
ties of states p a {oj) ~ Pa> such that 



G£(Mi)ffa (*i - t) ~ iG&&t)Pa, (26) 



with i = 1 for a — L and i = N for a = R. The 
expression (|14|l for the dc current leads to 

-\w R \ 2 p R (fi>)\pft{w)-JS N (u>)], (27) 

where we have defined the dc component of the gen- 
eralized densities of states at the sites where the time- 
dependent potentials act: 



i 

TO 



dt2Im[G? i (t,oj) 



(28) 



and the dc components of the spectral densities of occu- 
pation at those sites, pi( w ) : 

p» = - / dt{\G^(t,co)\ 2 \w L \ 2 p L (u;) 

TO Jo 

+\G? N (t,v)\ 2 \w R \ 2 p R ( U ), 



i r° 

m = -*— / dtGfifat) 
T o Jo 



du> 
2^ 



/Hft( W ),(29) 



being m the density of particles at the site i. In an equi- 
librium system, these two spectral functions coincide, i.e., 
Pi(uj) = pf c (u)) and are positive defined functions, but 
in a time dependent problem these two functions dif- 
fer in general. While ~p~i(to) is a positive defined func- 
tion, in a time-dependent problem pi (t, ui) may change 
sign as a function of u>. Actually, it is clear from Eq. 
(|27|l that the violation of the equivalence between these 
two spectral functions, is at the heart of the existence 
of a non- vanishing dc current. Another necessary con- 
dition is the breaking of left-right symmetry. This can 
be accomplished statically by, for example, considering 
Pl(w) 7^ p R (to) or wl ^ w R , but also dynamically, by 
recourse to a finite phase lag S ^ in the Hamiltonian 
(|25p. In what follows, we focus on the latter case. 



B. Results 

We consider the symmetric array of two barriers de- 
scribed in the previous subsection, which are placed be- 
tween two identical reservoirs with a large bandwidth 
Pl(w) — p R (u>) — P (u>) with the form defined in 



and \wl\ 2 — \w R \ 2 — \wq\ 2 . Due to the symmetry of the 
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FIG. 5: (Color online) The transmission function T(uj) for a 
pumping amplitude V — 0.2 a phase lag S = tt/2 and fre- 
quencies f2o = 0.01,0.3 (top and bottom panels). The el- 
lipse in the lower panel encloses a region where resonance is 
achieved. Other parameters are Eb = 1, jwo| 2 = 0.5, W = 4 
and N = 20. 



densities of occupation p 1 (to) and p N (oj). In this way, 
the transmission function l|27|l simplifies to 



(30) 



and it is found that the current behaves like J oc (njv — 
m). This is very suggestive, since it is natural to asso- 
ciate a difference in the local density of particles with a 
static potential drop. In the present case, we can imagine 
that the pumping with a phase lag induces an effective 
potential drop V e /f between the sites 1 and N, which 
causes a difference in the particle population at these 
two sites, and a current J oc V e ff. 

An important situation corresponds to the case of weak 
pumping amplitude V. A perturbative solution of the set 
© to the lowest order in V leads to 

^ e -<(*+ n o*)GS lil ( W + n )GS jn (a;) + 



(31) 



Using these expressions for the Green functions, making 
use of the fact that for the symmetric device G^^lj) = 
Gff N (u)) and G® n (uj) = G N1 (w) and keeping terms to 
the lowest non- vanishing order in V , we find 



TV) 



2V 2 \w 1 4 sm(6) [p a (u)] 2 Re{G° n (w) [G a 1N (u)}*} 
xIm{G° n (Lo + Sl )[G° 1N (u + n )}* 
-G° n (oj-n )[G 1N (co-n )]*}. (32) 



setup. 



fft{u>) 



pjy(uS), while 5 causes differences in the 



In order to have a better insight on the behavior of 
the function T(u>), we present some results in Fig. [5] for 
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FIG. 6: (Color online) The dc current J as a function of 
the chemical potential of the electrodes [i. Top and bottom 
panels corresponds to couplings to the macroscopic reservoirs 
|uio| 2 = 0.1 and |wo| 2 = 0.5, respectively. The pumping am- 
plitude is V = 0.2. Plots in thin solid black, dashed red, 
dot-dashed blue and thick solid magenta lines correspond to 
fio — 0.01,0.3,0.45,0.6, respectively. The spacing between 
two consecutive levels is referred to as A. Other parameters 
are as in Fig. 



a phase lag 5 — ir/2. In the figure, we show the exact 
T(uj), calculated from eq. I|17|) . As the bandwidth of the 
macroscopic leads (2W = 8) is large compared to the 
spectral width of the central piece (~ 4) , the wide-band 
description for the reservoirs applies and we have checked 
that T(u) is well reproduced by equation l|3U|l . The upper 
panel of Fig. [3] corresponds to small pumping amplitude 
V, small pumping frequency f2o and 6 = n/2. The trans- 
mission function shows a peak-antipeak pattern with a 
separation A that coincides with the energy distance be- 
tween two consecutive energy levels of the central struc- 
ture. The peak-antipeak behavior of T(uj) indicates that, 
as the chemical potential increases and covers a level, 
electrons with lower energies within the linewidth are al- 
lowed to travel from the left to the right while the ones 
with higher energy travel in the opposite direction. The 
corresponding behavior of the dc current as a function of 
/i is shown in thin black line in the upper panel of Fig. 
and it consists in a succession of small peaks suggesting 
that as the chemical potential is increased covering an 
energy level of the central system, a conduction channel 
is enabled and a net current flows between the two reser- 
voirs. The sign of the current is consistent with the one 
expected from intuitive considerations for this value of 
the phase lag on the basis of the following adiabatic pic- 
ture: the first barrier lowers its effective potential during 
the part of the cycle where the second potential grows, 
favoring the incoming of electrons from the left reservoir 
inside the well. In the other part of the cycle, the voltage 
of the second barrier gets lower, helping the electrons to 
tunnel from the well towards the right reservoir. 

An interesting situation takes place when the fre- 



quency is resonant, i.e. £lo = A. In this case, two 
neighboring electronic levels of the central device are ex- 
pected to be mixed by the pumping potentials. The re- 
gion marked in the lower panel of Fig. [S] satisfies the 
resonant condition and the function T(u>) consists of a se- 
quence of antipeaks. The corresponding current is shown 
in dashed red lines in the upper panel of Fig. and ex- 
hibits plateaus within this region. The sign corresponds 
to a net electronic flow from the right to the left suggest- 
ing that the quantum interference due to the mixing of 
levels causes a sign reversal in comparison to the situa- 
tion observed for pumping frequencies Qq < A. Also note 
in Fig. El that when the resonant condition is £lo = 2A, 
the current recovers the sign of the situation Slo < A. 
For weak pumping amplitude V, a further examination 
on the origin of this interference is possible on the basis 
of the approximate solution given in Appendix B. The 
analysis presented there indicates that, at resonance, the 
effective phase lag between the two potentials is S + jir, 
being j — 1 the number of energy levels of the structure 
between the two interfering ones. On the other hand, due 
to the symmetry of the problem, J(6 = (— I) 3 J(S), 

which complements the argument to explain why a shift 
in the phase lag may cause a change in the sign of the 
net current. The lower panel of Fig. gives an idea 
of the effect of the changes in the linewidth of the cen- 
tral piece introduced through an stronger coupling to the 
reservoirs. As expected, these effects are more important 
at resonant frequencies since they help in the mixing of 
the electronic levels. For high energies, close to the up- 
per edge of the spectrum (uj ~ 2), the energy levels are 
closer (A is smaller). Therefore, the resonant frequencies 
strongly differ from the ones at lower energies and sign 
reversals of J are observed as a function of /i. 

The behavior of J as a function of the phase lag is 
shown in Fig. [7] for a selected value of the chemical po- 
tential [i. All the plots of the upper panel correspond to 
a small pumping amplitude and can be fitted by a func- 
tion oc sin(<5), in full agreement with eq. I|32l) . The lower 
panel corresponds to a higher V and deviations from this 
behavior are found. For example, the plot in dashed lines 
corresponds to a function of the form J = Asin(<5) and 
fits the behavior of J for f^o = 0.6 only in a neighbor- 
hood of S = 7T. Similar deviations are observed for other 
pumping frequencies. 

The behavior of J as a function of the coupling to the 
reservoirs \wo\ 2 is shown in Fig. [5] for fixed /i, with dif- 
ferent frequencies and pumping amplitudes. The higher 
panel corresponds to a low pumping amplitude and it is 
drawn in log-log scale in order to observe details for very 
low coupling to the reservoirs. In most of the cases, the 
current vanishes as the coupling to the reservoirs tends 
to zero following a law J oc |wo| 4 a s suggested by 
The corresponding reference slope is indicated in the fig- 
ure for comparison. In some cases, slight deviations from 
this law are observed. One example is the plot in red 
dashed lines shown in the upper panel which corresponds 
to a resonant frequency. The origin for such departures 
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FIG. 7: (Color online) The dc current J as a function of 
the phase lag 5 for fi = and |wo| 2 = 0.5. Upper and lower 
panels correspond to pumping amplitudes V = 0.2 and V = 1, 
respectively. Circles, squares, triangles, diamonds and stars 
correspond to Slo = 0.01, 0.05, 0.3, 0.45, 0.6, respectively. The 
function J = 0.0185 sin (5) is indicated in dashed lines in the 
lower panel. Other parameters are as in Fig. |S] 



should be found in the fact that the functions G°- (w) tend 
to be singular as the coupling to the reservoirs tends to 
vanish. The case of a larger pumping amplitude is il- 
lustrated in the lower panel of Fig. [8] A notable feature 
observed in the latter case is that changes in the strength 
of the coupling to the reservoirs may introduce changes 
in the sign of J, like in the case of the plot corresponding 
to f2o = 0.75. This is because the coupling to the reser- 
voirs contributes to enhance the quantum interference 
between energy levels. For very small |u>o| 2 , the behavior 
is similar to the one observed for weak V, illustrated in 
the higher panel. 

Figure [5] shows J as a function of the pumping fre- 
quency f2o for different values of the parameter |u>o| 2 . 
The first feature to note is the structure of minima and 
maxima corresponding to resonant frequencies causing 
interference between nearest and next-nearest neighbor 
energy levels. The change of sign between resonances is 
consistent with the arguments of Appendix B and the 
discussion related to Figs. El and [fj] Another issue worth 
mentioning is the linear behavior at very small pumping 
frequencies, as can be inferred from an expansion of the 
low V transmission function IMl'li in powers of Qq. 

To finalize, we show in Fig. J as a function of the 
pumping amplitude V. For chemical potentials bellow 
the spectral edge, the behavior is consistent with J oc V 2 
within a rather wide range of V. This can be better 
observed in the log-log plot of the inset. This law is in 
agreement with the low V transmission function (|32|l . 
Note that this behavior is followed at resonance within a 
wider range of V (see the plots in blue dot-dashed and 
thick magenta lines in the left panel). The right panel 
shows the situation for chemical potentials close to the 



FIG. 8: (Color online) The dc current J as a function of the 
contact parameter |uio| 2 for 6 = it/2. The upper panel is plot- 
ted in log-log scale and corresponds to a pumping amplitude 
V — 0.2 and chemical potential fi = 0. The reference slope 
corresponding to J oc |u>o| 4 is plotted in light blue dotted line. 
The lower panel corresponds to (i = 1 and pumping ampli- 
tude V = 1. Plots in thin solid black, dashed red, dot-dashed 
blue and thick solid magenta and thick dashed-dot orange 
lines correspond to Qq = 0.1,0.3,0.45,0.6,0.75, respectively. 
Other parameters are as in Fig. 



0,002 
0,0015 

0,001 
0,0005 


-0,0005 

-o.ooi 

-0,0015 



1 1 


' 1 1 1 


' 1 1 


1 





I 1 



0,03 
0.02 
0,01 


-0,01 
-0,02 



1 


1 ' 


I'll 

S - 








1 


1,1,1 


1,1 







0.1 



0.3 



FIG. 9: (Color online) The direct current J as a function of 
the pumping frequency Qo for S = n/2, chemical potential 
fj, = 0, and contact parameter |u>o| 2 = 0.1,0.5, 1 (black solid, 
red dashed and blue dot-dashed lines). Upper and lower pan- 
els correspond to pumping amplitudes V = 0.2 and V = 1, 
respectively. Other parameters are as in Fig. 



spectral edge. In this case, the behavior for very small V 
remains in agreement with the law predicted by i|32f) . For 
some pumping frequencies, however, J strongly deviates 
from this behavior at moderate V, in some cases showing 
sign reversals. Such dramatic chan ges h ave also been ob- 
served in molecular driven ratchets I14LI15I. Their source 
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FIG. 10: (Color online) The direct current J as a func- 
tion of the pumping amplitude V for S — tv/2 and fio = 
0.01, 0.1, 0.3, 0.6 (black thin, red dashed, blue dot-dashed and 
magenta thick lines, respectively) and |?a>o| 2 = 0.5. Left and 
right panels correspond to y, = and n = 2, respectively. 
The inset in the left panel contains the same data of the main 
panel in log-log scale. The light blue dashed line indicates the 
slope for J cx V 2 . Other parameters are as in Fig. 

must be found in the higher density of electronic levels at 
higher energy and the possibility of interference between 
more than two electronic levels as the pumping ampli- 
tude increases. As in the case of only one harmonically 
time-dependent potential, the behavior of T(uj) within 
the high energy region has a much higher sensitivity to 
the particular values of parameters. Hence, the predic- 
tion of the sign of the current when the chemical potential 
tunes this spectral region turns much more difficult. 



V. SUMMARY AND CONCLUSIONS 

We have presented a general approach based on non- 
equilibrium Green functions to study transport phenom- 
ena originated in time-periodic potentials applied on 
quantum systems described by tight-binding Hamiltoni- 
ans without many-body interactions. The present treat- 
ment allows for the exact solution in problems with sev- 
eral time-dependent potentials provided that the oscil- 
lating frequencies are commensurate (i.e. a multiple of 
an elementary frequency Qq). In addition, the present 
treatment is valid for arbitrary amplitudes and oscillat- 
ing frequencies of the time-dependent potentials. 

We have employed the general treatment of section 
II to study two simple problems of quantum pumping 
in a one-dimensional model for a double barrier struc- 
ture connected to left and right reservoirs. We have first 
considered only one pump acting on one of the barri- 
ers. This case is interesting because an explicit analyt- 
ical expression for the retarded Green function can be 
found. We have shown that the existence of a net elec- 
tronic transport depends exclusively on the structure of 
the environment of the pumping center and two simple 



conditions must be fulfilled: the geometrical arrangement 
must not have spacial inversion symmetry with respect 
to the pumping center and it must have resonant levels. 

The second case we have considered corresponds to 
pumping potentials oscillating with the same frequency 
and a phase-lag between them, acting at both of the 
barriers. This kind of operational arrangement of ac- 
potentials is just the one used in the experiment of Ref. 
|2|. In particular, we have shown that for weak pump- 
ing amplitude, the net current behaves like J oc sin(5) 
as a function of the phase lag, in agreement with the ex- 
perimental work. Another interesting feature is that for 
reservoirs with wide bands (as is often the case in double 
barrier structures in semiconductor junctions), the cur- 
rent can be related to the difference in the charge density 
at the pumping points. 

We have stressed the fact that all electrons below the 
Fermi energy of the reservoirs contribute to the current 
generated by pumping mechanisms. We have also de- 
voted some effort to understand the conditions that de- 
termine the direction of the net current. We have found 
that in the two cases analyzed, this direction coincides 
with the one predicted by a "naive" adiabatic picture in 
the situations where quantum interference does not play 
a role. This condition seems to be more easily achieved 
in the case of only one pumping potential. Instead, for 
the case of two pumping potentials it is achieved for not 
too strong pumping amplitudes, when the pumping fre- 
quency is smaller than the separation between consecu- 
tive energy levels. For resonant frequencies, sign reversals 
take place due to quantum interference of different elec- 
tronic levels. In both examples, there are regions close to 
the spectral edge where the transmission function exper- 
iments many changes of signs. In this region, the details 
strongly depend on the parameters. We have also shown 
that the current behaves linearly in f2o when the pump- 
ing frequency is small and proportional to V 2 for low 
pumping amplitudes. At resonance, the latter behavior 
extends to a range of V beyond the one expected from 
arguments based in perturbation theory. 

The investigation of the pumping effect in systems 
with annular topology and the combination of the time- 
dependent effects with electron-electron, electron-phonon 
interactions and disorder is left to the future. 
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APPENDIX A: EVALUATION OF THE 
RETARDED GREEN FUNCTION FOR A SINGLE 
HARMONICALLY TIME-DEPENDENT 
POTENTIAL. 



APPENDIX B: APPROXIMATE SOLUTION FOR 
THE RETARDED GREEN FUNCTION FOR TWO 
HARMONIC POTENTIALS AT RESONANCE. 



We present the solution of the set (JSJ for the case of 
only one periodic potential. In this case, it is possible 
to obtain an analytical expression for Gf"^, u>) through 
a recursive procedure. After some algebra, it is found 
an expression with the structure ©, with the following 
coefficients: 

&,!(*, w) - (y) fe g e// M J] g { - m) {uj - mfio), k > 



(Q,w) 
where 



m—1 
-k 



(Al) 



E e// (w) 



-£l - S e //(w) 



E°H + (y) L9 (1) (^ + «o) 

+ ff (-i)( w _n )]. 



(A2) 



The "bare" self-energy E°(w) = |w L | 2 .gf' ( w ) + l w fl| 2 5il( w ) 
represents the environment and is completely defined 
from the free Green functions of the pieces at the left 
(L) and at the right (i?) of the pumping center. The 
function g^ m \uj + mfl ) can be expressed as a continued 
fraction defined from the recursion relation 

[g {m) {uo + mflo)}- 1 =lu + mflo - £i - + mQ ) 



-(-) g( m±1 Hcj + ( m ±l)n ), 



(A3) 



where +, — corresponds to to > and to < 0, respec- 
tively. In practice, a cut-off is introduced such that 
g (±K \uj±Kno) = [u±Kfl -FP(LJ±Kn )]- 1 , being K 
large enough in order to satisfy that |iff2o| is much larger 
than the absolute value of the energy for which the bare 
Green function [u> — S (tj)]^ 1 has non-vanishing spectral 
weight. 

Using the fact that 

Im[£ M (0,^)] = -Im[£ "(w)]|$ e// M| 2 , (A4) 
and the definition of E 6 *-' (ui), it is found 

lm[ai,i(0,w)] = -^Im[Z°(uj - kSl )]\g hl (k,u)\ 2 . 
k 

(A5) 



Let us start by noting that the retarded Green function 
for a tight-binding chain with hopping w and length N 
with open boundary conditions at both ends is: 



k 



sin(fcZ) sin(fcm) 
lu - e k + ir\ 



(Bl) 



being r\ = 0+, e k 
n = l,...,N. 



-w cos(fc), and k — nil/ (N + 1), with 



We assume that the resonant energies e k and phases 
of the wave function of the double barrier structure con- 
nected to the reservoirs are approximately those of an 
open tight-binding chain of the same length and we pro- 
pose the following ansatz for the retarded Green function 
evaluated at a resonant frequency e k 



G °i m ( e k) ~ 7fc sin(fcl) sin(fcm). 



(B2) 



We have verified that this ansatz reproduces with rea- 
sonable accuracy the phases of the Green function in the 
present model when the coupling to the reservoirs is not 
too weak. 

Let us consider weak pumping amplitude and a pump- 
ing frequency such that 0,$ = e kl — e k2 , being fci — k% — 
jn/(N + 1) being j a positive integer. The set <|31ll re- 
duces to 



G* n (*,e fcl ) ~ G^,„(e* 1 ) + ^e i (' +no *^ il ( ejfca )GS jn ( efcl ) 
+ _ e 4f2 o t G^j V (e fc2 )G^ n (e fcl ), 

+ je- m ° t G° mfN (e kl )G%, n (e k2 ). (B3) 



The replacement of the ansatz <|B2[) in the above equa- 
tions reveals that there is a phase equal to jn between 
the terms proportional to y e ± l ( s + n at) - m re i a tion to the 
ones proportional to Ve ±lQ ~ ot in the above expressions for 
G^ n (t,e ki ). This indicates that the effective phase-lag 
between the two potentials is S + jn. 
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